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In a classical ideal plasma, a magnetic field is known to reduce the heat conductivity perpendicular 
to the field whereas it does not alter the one along the field. Here we show that, in strongly correlated 
plasmas that are observed at high pressure or/and low temperature, a magnetic held reduces the 
perpendicular heat transport much less and even enhances the parallel transport. These surprising 
observations are explained by the competition of kinetic, potential and collisional contributions to 
the heat conductivity. Our results are based on first principle molecular dynamics simulations of a 
one-component plasma. 
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Strongly coupled plasmas have become a focal point of 
plasma research in recent years as several experimental 
setups such as dusty plasmas [1, 2], trapped ions [3], and 
ultracold plasmas [4] have become available to study this 
unusual state of ionized matter in which the potential 
energy exceeds the kinetic energy. Similar plasma condi¬ 
tions are thought to exist, e.g., in white dwarf stars or in 
the outer layers of neutron stars [5, 6]. 

Often, these plasmas are subject to strong magnetic 
fields, as is the case in neutron stars, in magnetized tar¬ 
get fusion experiments [7], or in recent dusty plasmas 
experiments [8]. 

The interplay between strong particle interaction and 
a strong magnetic field gives rise to a rich dynamics [9], 
including complex wave spectra [10-12], generation of 
higher harmonics [13, 14], and long-lived metastable 
states [15]. The diffusive particle transport has been an¬ 
alyzed in two and three dimensions [16-18] and in binary 
plasmas [19] with the conclusion that the self-diffusion 
coefficient is strongly reduced by both the strong cou¬ 
pling and the magnetic field. 

Of key importance for plasmas and their applications 
is their capability to conduct and transport heat. This 
plays a crucial role, e.g., in the plasma confinement in 
modern fusion scenarios [7] or for the cooling rate of 
neutron stars [20]. Heat transport is well understood in 
magnetized weakly coupled plasmas since the pioneering 
work of Braginskii [21], and was studied in strongly cou¬ 
pled unnragnetized systems, e.g., in dusty plasma exper¬ 
iments [22-24], by kinetic theory [25], and by computer 
simulations [26-31]. However, until now a systematic 
study of heat transport in magnetized strongly coupled 
plasmas has not been carried out [32]. It is the goal of 
this work to fill this gap. By performing first principles 
molecular dynamics (MD) simulations of one-component 
plasmas (OCP) we show that the heat conductivity A 
radically differs from the one in weakly correlated high- 
temperature plasmas [21, 34]: a) the reduction of Aj_ (the 
component perpendicular to B) with increasing magnetic 


field strength B is much slower and Ax approaches a non¬ 
zero asymptotic value, and b) the parallel component A|| 
increases with B instead of remaining constant. 

Theory and simulation approach. In a magnetized 
system, the energy flux j is related to the tempera¬ 
ture gradient via the thermal conductivity tensor as 
j a = —\ a p (VT)^, where the latter has three indepen¬ 
dent components (we assume B || e z ), 


/ A± A x 0 \ 

A = -A x A ± 0 . (1) 
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The components A|| and Ax describe field-parallel and 
cross-field heat transport, respectively, and converge to 
the scalar heat conductivity A = Ao as B —> 0. The 
off-diagonal term A x is the analog of the Hall effect (the 
so-called Riglri-Leduc effect [35]) and vanishes for B —► 0. 

A microscopic approach to X a p is provided by the 
Irving-Kirkwood expression [36, 37] for the heat flux, 
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which consists of kinetic, potential, and collision contri¬ 
butions (first, second, and last term). 

The Green-Kubo formula then relates A to the auto¬ 
correlation function [38], 


A q /3 = lim 1 [ (j a (t)j 0 (O))dt, (3) 
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where the integral is over the energy autocorrelation 
function. According to Eq. (2), each component of \ a p is 
the sum of three direct and three cross-correlation terms. 

To compute \ a p for a typical strongly coupled OCP, 
we use a screened Coulomb (Yukawa) potential = 
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FIG. 1. Heat conductivity in an unmagnetized Yukawa sys¬ 
tem at k = 2 as a function of F. Shown are the total conduc¬ 
tivity and its kinetic, potential, and collisional parts, “cross” 
denotes the sum of the three cross-correlation terms. 


Q 2 /r x exp (— kr) for the pair interaction with the in¬ 
verse screening length k and apply a homogeneous mag¬ 
netic field. The system is then completely character¬ 
ized by three dimensionless parameters: the coupling 
parameter T = Q 2 / (ksTa), the screening parameter 
k = ak , and the normalized magnetic field strength 
/ 3 = uj c /u}p oc B , where a = [3/ (Ann)] 1 ^ 3 is the Wigner- 
Seitz radius, u> p = [47t Q 2 n/m} 1 ^ 2 is the plasma frequency, 
and lo c = qB/(mc ) is the cyclotron frequency (Q,T, and 
c are the charge, temperature, and speed of light, respec¬ 
tively). For later use we also define a = t co iui p , where 
t co i is the collision time. 

We performed extensive MD simulations for magnetic 
fields [39, 40], using N = 8192 particles in a cubic 
box with periodic boundary conditions and a value of 
n = 2 which is representative for a strongly coupled OCP. 
Data collection begins after an equilibration period and 
takes place under microcanonic conditions for a time of 
tojp = 2.5 x 10 5 . Since the evaluation of the integral (3) 
is limited to a finite time r in practice, its value is deter¬ 
mined by averaging Eq. (3) over uj p t £ [1000,2000] [41]. 
For each data point, 50 separate simulations with differ¬ 
ent initial conditions are averaged for a total measure¬ 
ment time of ui p t = 1.25 x 10'. We report the standard 
error of the mean of these simulations; unless shown, the 
error bar for all values is (much) smaller than the symbol 
size. Data for the heat conductivity are given in units of 
nkBU p a 2 . 

Unmagnetized system —We start with an isotropic, un¬ 
magnetized system where \ a p is diagonal. Figure 1 re¬ 
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FIG. 2. Elements of the heat conductivity tensor at different 
magnetic field strengths as a function of F (top), and the value 
at /? = 1 relative to the held-free value Ao (bottom). 


veals an interesting non-monotonic dependence of A on F, 
which is one of the hallmarks of strongly coupled systems. 
The present results represent a significant improvement 
over previously available data and facilitate a more de¬ 
tailed analysis. 

In particular, the decomposition (2) uncovers the ori¬ 
gin of this dependency: The kinetic contribution to A 
corresponds to the (material) transport of kinetic energy 
associated with the movement of a particle between col¬ 
lisions. It falls off, approximately as F -3 / 2 , since both 
the mean free path and the thermal energy per particle 
decrease with T [42]. Likewise, the potential contribution 
is associated with the material transport of potential en¬ 
ergy and increases with T as the average potential energy 
grows. In the region of T = 30... 200 this growth slows 
down, since the mean free path decreases with T, and the 
system dynamics transform into caged motion. Finally, 
the collisional part increases with T up to the crystalliza¬ 
tion point (T c « 420) where heat transport is dominated 
by phonons [43]. 

Cross-field heat transport, Aj_- In the presence of a 
magnetic field, the T-dependence of A changes drasti¬ 
cally, as is shown in Fig. 2. Consider first the field- 
perpendicular contribution, Aj_. As in a weakly cou¬ 
pled plasma, A_l is reduced by the magnetic field, how¬ 
ever, with increasing F this reduction becomes less pro¬ 
nounced. This is readily explained by the change of the 
governing heat transport mechanisms in different cou¬ 
pling regimes: For small T, material transport of kinetic 
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FIG. 3. Contributions to the heat conductivities Ax, Ay, and 
A x at /3 = 1. “cross” denotes the sum of the three cross¬ 
correlation terms. 

energy dominates (see Fig. 1) which is greatly reduced 
perpendicular to B by the cyclotron motion of the parti¬ 
cles. On the other hand, at large F, the thermal conduc¬ 
tivity is mainly due to collisions between particles (Fig. 1) 
whose frequency and effectiveness are only weakly re¬ 
duced by the magnetic field. Figure 3a shows how these 
different effects give rise to the observed T-dependence of 
Aj_ at P = 1. 

Let us now analyze the dependence of Ax on the mag¬ 
netic field strength. For weakly coupled systems, both 
classical transport theory [21] and a hydrodynamic anal¬ 
ysis [44, 45] predict a decay of Aj_ scaling as Ax ~ /3 -2 , 
whereas the parallel conductivity is independent of the 
field strength, Ay ~ j3°. The corresponding simulation re¬ 
sults for strongly correlated systems are shown in Fig. 4 
for the cases T = 5 and F = 40. 

As observed before, Aj_ decreases with /?; however, it 
does not vanish, but approaches a finite value. The anal¬ 
ysis shows that this is due to a residual heat transfer 
via collisions, see Fig. 4: even if particles are unable to 
move perpendicular to B, heat is still transferred by col¬ 
lective modes of the plasma, in particular, the ordinary 
shear mode and the upper and lower hybrid modes [11], 
a mechanism similar to heat transfer via phonons in crys¬ 
tals. 

Thermal Hall effect A x - The off-diagonal tensor com¬ 
ponent (cf. Figs. 2, 3b, and 4) is a special case as it 
emerges only in the magnetic field. However, an increase 
of B leads to smaller Larmor radii, which decreases the 
efficiency of this transport mode. These competing ef¬ 
fects result in a non-monotonic behavior with a maxi¬ 
mum around = 0.1... 0.5, cf. Fig. 4. With increas¬ 
ing coupling A x decreases but approaches a finite value, 
cf. Fig. 3b, which is explained by the decrease of the Lar¬ 
mor radius with T (decrease of temperature) [46]. 

Field-parallel transport , Am—F igure 2 shows a striking 
result: heat conduction parallel to B is enhanced by the 
field. This result is surprising, since it is in contrast to 
the behavior of weakly coupled plasmas. But even for 



FIG. 4. Field dependence of the heat conductivities. Thick 
lines represent the total values of Ax, Ay, and A x and thin 
lines of the same color show the respective kinetic, potential 
and collisional contributions. 

strongly coupled plasmas this is an unexpected field ef¬ 
fect as it contradicts the behavior of the diffusion coef¬ 
ficient parallel to B which has been found to decrease 
monotonically with B [17]. 

To clarify the origins of this enhancement, consider 
again the different contributions to Ay, Fig. 3c. There is 
a drastic increase of the collisional contribution which is 
easy to understand: collisions parallel to B are becoming 
more effective with increasing field strength, since parti¬ 
cles are unable to avoid one another in lateral direction 
after the initial approach. This leads to a growing inter¬ 
action time and, subsequently, an increased energy ex¬ 
change. This is particularly important for moderate cou¬ 
pling (r < 10), where infrequent binary collisions domi¬ 
nate, and less so in highly coupled systems, where caged 
motion already results in highly efficient collisions. 

The potential contribution to Ay is likewise enhanced 
by the field. Here, a competition between two processes is 
observed: On the one hand, the mobility of the particles 
(i.e., the diffusion coefficient) along B is reduced by the 
field [17]. On the other hand, the reduced interaction 
with neighboring particles in the cross-field plane enables 
the particles to retain energy for a longer time, resulting 
in a net increase of the heat conductivity. 

Similar considerations apply to the kinetic contribu¬ 
tion. 

The interplay of these processes, together with the 
varying relative importance of the kinetic, potential, and 
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FIG. 5. Relative heat conductivities (top) and effective total 
heat conductivity A cfi = TrA (bottom) compared to classical 
transport theory for weak coupling, gray curves [34]. To fa¬ 
cilitate the comparison, values of a = 2.5 (a,c) and a = 1.24 
(b) were used (see text). 


collisional contributions, work to increase A|| over Ao for 
the whole range of coupling strengths considered. Note 
that the enhancement of Ay is weaker for smaller T, and 
Ay —> Ao for F < 1 (see Fig. 2), confirming consistency of 
our simulations with the weak coupling limit. 

Consider now the dependence Ay(/3). As Fig. 4 shows, 
A|| approaches a maximum value as /? increases. This 
occurs regardless of whether transport is dominated by 
material transport (F = 5) or by collisions (r = 40). 
The physical reasons for this is that both the increase of 
energy retention and of the efficiency of collisions, which 
are the driving mechanisms behind the growth in Ay, have 
upper limits, i.e., complete energy retention and complete 
collisional energy transfer. 

Test of classical transport theory —Our simulation re¬ 
sults allow us to determine the applicability limits of the 
weak coupling theory of Braginskii [21, 34] which pre¬ 
dicts the scalings A_l(/3) ~ a~ 2 /3~ 2 and A x (/3) ~ a -1 /? -1 
for large ft. A comparison with the strong coupling data 
at hand shows semi-quantitative agreement for A_l and 
A x for moderately coupled plasmas, T < 5, when a(r) is 
regarded as a free parameter. Figure 5 shows the best fit 
for r = 5 (a = 2.5) and T = 40 (a = 1.24) [47] which in¬ 
dicates excellent (poor) agreement in the former (latter) 
case. At the same time, this extension of weak coupling 
theory neither captures the finite asymptotics of Ay for 
large /3, nor can it reproduce the increase of Ay with /?. 

Summary —We have presented new high-precision data 
for the thermal conductivity of unmagnetized and mag¬ 


netized Yukawa one-component plasmas. The decompo¬ 
sition of A into the different modes of transport (kinetic, 
potential, and collisional) has enabled us to elucidate the 
origin of the well-known non-monotonic dependence of 
A(r) at zero magnetic field for the first time. 

At finite magnetic fields, the cross-field thermal trans¬ 
port is reduced, as in weakly coupled plasmas; however, 
the decay is close to 1/B at intermediate field strengths 
and approaches a finite value at large fields which is due 
to the collisional mode. On the other hand and, con¬ 
trary to what one might expect by extrapolating from 
the weakly coupled case, the field-parallel heat conduc¬ 
tivity is enhanced by the field in strongly coupled plas¬ 
mas. This is due either to an enhanced energy retention 
in systems where kinetic transport dominates or due to 
a modification of the collision process in systems domi¬ 
nated by this transport mode. 

Finally, to understand the implications of our re¬ 
sults for strongly coupled plasma applications we con¬ 
sider a heated plasma cube of volume V. The energy 
loss Jt(T,/ 3) through the surface (and, thus, the tem¬ 
perature change) will be proportional to [cf. Eq. (1)] 
A eff = TrA = Ay + 2Ay. At weak coupling, Jy can be re¬ 
duced by increasing the magnetic field strength—a stan¬ 
dard tool in many plasma applications—to one third of 
the field-free value (only Ay remains), see Fig. 5c. How¬ 
ever, in strongly coupled plasmas, our results indicate 
that using a magnetic field to reduce heat losses is much 
less effective: for large coupling (cf. curve for T = 40) 
there is almost no reduction whereas, for moderate cou¬ 
pling, (cf. curve T = 5) losses will be reduced for weak 
fields, but they will be substantially enhanced in strong 
fields. Our results have important consequences for the 
heat transport of strongly coupled magnetized plasmas. 
One example are the cooling rates of neutron stars and 
magnetars. The second concerns inertial confinement fu¬ 
sion setups that may, in fact, benefit from the reduction 
of the heat conductivity and diffusion [17] in strongly 
coupled plasmas, cf. Fig. 1. In this case, further sub¬ 
stantial reduction of A by a magnetic field is possible but 
requires the use of asymmetric fuel geometries with a low 
surface area in field direction. 
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